Rho Function

private function Rho(lat, a, e) result(radius)

Compute the radius of curvature of the ellipsoid in the plane of the meridian at given latitude

Arguments

Type IntentOptional Attributes Name
real(kind=float), intent(in) :: lat

latitude [radians]
real(kind=float), intent(in) :: a

semimajor axis [m]
real(kind=float), intent(in) :: e

exentricity

Return Value real(kind=float)


Source Code

FUNCTION Rho &
!
(lat, a, e) &
!
RESULT (radius)

IMPLICIT NONE

!Arguments with intent(in):
REAL (KIND = float), INTENT (IN) :: lat !!latitude [radians]
REAL (KIND = float), INTENT (IN) :: a !!semimajor axis [m]
REAL (KIND = float), INTENT (IN) :: e !!exentricity

!Local variables:
REAL (KIND = float) :: radius

!------------end of declaration------------------------------------------------

radius = a * (1. - e**2.) / (1. - e**2. * SIN(lat)**2) ** 1.5

END FUNCTION Rho